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Abstract 

A novel formal equivalence between thermal averages of coherent properties (e.g. conductance), 
and time averages of a single wavepacket arises for Fermi gasses and certain geometries. In the case 
of one open channel in a quantum point contact (QPC), only one wavepacket history, with wavepacket 
width equal to thermal length, completely determines thermally averaged conductance. The formal 
equivalence moreover allows very simple physical interpretations of interference features surviving 
under thermal averaging. Simply put, pieces of thermal wavepacket returning to the QPC along inde- 
pendent paths must arrive at the same time in order to interfere. Remarkably, one immediate result of 
this approach is that higher temperature leads to narrower wavepackets and therefore better resolution 
of events in the time domain. In effect, experiments at 4.2 K are performing time gated experiments 
at better than a gigahertz. Experiments involving thermally averaged ballistic conductance in 2DEGS 
are presented as an application of this picture. 

1 Introduction 

Quantum coherence and control has become a major theme in device physics, with applications in quan- 
tum information and computing as possible rewards for success. The quantum world offers new phenom- 
ena involving interference and entanglement which are not incorporated into present day devices. The 
ability to store, process, and retrieve information coherently may lead to qualitatively new technologies. 

Finite temperatures tend to desttoy quantum effects, especially if kT is much larger than the mean 
level spacing of the system. However thermal averaging is not decoherence; it is not the same as entan- 
glement of the "system" with many other degrees of freedom. For, if we thermally average over states 
which are individually coherent with respect to some degrees of freedom of interest, that coherence re- 
mains in some sense. Well known examples from optics include the "white light interferometer", which 
operates even with blackbody radiation to give high contrast interference fringes. 
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Coherent wavefunction descriptions of matter waves give way to incoherent density matrix descrip- 
tions at finite temperature. The common wisdom is that a single coherent wavefunction will not suffice 
when describing a system at thermal equilibrium. Although this is usually true, there is an exception 
which exists in principle for non-degenerate systems. Consider an operator A and a system with Hamil- 
tonian H with a non-degenerate set of eigenf unctions. Suppose a diagonal density matrix correctly 
describes the physical situation: 

/'T = ^Pn |V'n)(V'n I (1) 
n 

thermal average of the operator A is then {A)t = Tt:[ptA]. Now construct a coherent quantum wave- 
function 

H^) = XI exp[ixn] IV'n) (2) 

n 

Then the time average of A over the history of the wavepacket is then 

T 

<A>= lim - [ dt = y^pn{iPn\A\^n) = Tt[ptA] = {A)t, (3) 



i.e. the time average of the operator A over the "thermal wavepacket" cl){t) is the correct density matrix 
trace of A. The phases Xn may be chosen arbitrarily, or chosen so as to control the initial location of the 
wavepacket. Thus the history of a single, coherent wavepacket determines thermal average. 

There are several reasons why this is normally not very useful. First, it is usually not possible to 
construct a priori a wavepacket that has the desired pn's. For example, in a closed billiard any compact, 
analytically simple wavepacket will have wildly fluctuating Pn&- Second, very long time propagations 
may be needed to converge the time average, a costly and usually non-intuitive process. Third, time 
averaging won't resolve any exact degeneracy and cannot provide a proper incoherent average over the 
degenerate states. Since degeneracies abound in real systems, these systems are excluded (this would 
seem to exclude all scattering systems except possibly ID wires). Fourth, even if the p^s can be realized 
in a pre-determined wavepacket, their form is often such that the wavepacket is not intuitively or com- 
putationally useful. For example, neither the Boltzmann distribution nor the Fermi distribution gives a 
compact wavepacket if ksT is greater than the lowest characteristic energies of the system. 

Thus it is remarkable that all four of these problems moderate or disappear altogether for at least one 
very important experimental situation: determination of conductance G of electrons at finite temperature 
through quantum dots and devices coupled to thermal reservoirs via quantum point contacts (QPC's). 
Indeed for a single mode contact we shall show that a nearly Gaussian wavepacket of width equal to 
thermal length £t = yp/T^kBT where vp is the Fermi velocity, can be run to give thermally averaged 
conductance and other thermally averaged properties. 

Thermal averaging is normally thought of as diluting or washing out specific quantum information. 
Here, we show thermal averaging permits the deconstruction of the ordering of events in an ersatz time 
dependent experiment wherein a narrow electron wavepacket is injected into the device. The width of 
the ersatz wavepacket decreases with increasing temperature and the implicit time resolution improves. 
This raises the hope that finite temperatures might be harnessed to enhance quantum coherent control. 
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2 Thermal wavepackets 



Some time ago it was demonstrated that narrow waveguide channels exhibit quantized condutance, in 
accordance with the number of energetically allowed transverse modes of the guide QEl. Such conduc- 
tance quantized QPC's are now standard equipment in the field of two degree of freedom electron gas 
(2DEG) physics. We will mostly consider the limit in which or 1 channels are open (below or on the 
first conductance plataeu). Consider electrons at fixed energy E, using the language of scattering theory 
and the scattering matrix. Let t be the part of the scattering matrix describing transmission from left to 
right across the contact. Then, if the lowest transverse mode dominates the transmission (true from the 
tunneling limit and onto the first conductance plateau) the matrix T(£') = t{E)t{E)^ has dimension 
A^, where N is the number of modes in the left lead. However T{E) has matrix rank one, and has 
one nonzero eigenvalue ri (E) ; the remainder, — 1 of them, are 0. The channel which conducts with 
probability ti (E) is a particular- linear combination of the usual transverse modes of the leads on the 
left. This linear combination is given by the eigenvector corresponding to the nonzero eigenvalue ti (E) 
of T(£'). When T{E) is rank one, an arbitrary transverse mode (or incoming plane wave if we choose) 
impinging from the left produces the same state on the right as any other incident mode, apart from an 
overall complex amplitude. 

Consider a clean, adiabatic QPC with one open mode and therefore conductance 2e^ /h lying between 
two perfect reservoirs. We introduce various "imperfections" to the right of the QPC: scattering impu- 
rities, moveable obstructions (such as the potential hill induced by a charged scanning probe tipll3ll4l). 
charged gates making confined dot regions, donor potential fluctuations, etc. These structures typically 
have the effect of reducing the overall conductance below 2e^//i, because they cause some scattering 
back through the QPC. (However, if they induce resonances on the right side, this can increase the flow 
through the QPC, so that the conductance for example reaches 2e^ /h even in the nearly pinched off 
regime. See refs. ||5l|6l.) 

We wish to determine thermally averaged conductance, including the structures on the right . The 
conductance is given at a single energy by QUI 



For an adiabatic, reflectionless contact between two clean, open reservoirs, on the first conductance 
plateau, Ge = 2e^//i; i.e. ti{E) = 1 in that energy region. The thermally averaged conductance is 
given byjS] 



where fx is the Fermi distribution function at temperature T. Assuming that the second (excited) and 
higher transverse modes in the contact are energetically out of reach, the thermally averaged conductance 
Eq.|5lis an average over channel 1 only: 



We denote by ipi{x,y,E) the total scattering eigenstate with the incoming boundary condition of 
precisely the eigenchannel of t which propagates with with probability ri (E) through the QPC (and 



Ge = -i-Tr[t{E)m^ = —n{E) 



(4) 




(5) 




(6) 
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which gives the 2e^//i ti{E) conductance, before the structures on the right were introduced). All the 
other degenerate channels orthogonal to ipi{x,y, E) have zero transmission. Consider the wavepacket 

^PT{x,y,Q) = j dEariE) e'^(^) M^.y.E). (7) 



where ariE) = hr^ ^ —m {dfT{E)/dE)/2'K, which gives a unit normalized wavepacket. The phase 
ip{E) can be chosen so that t/Jrix, y, 0) represents a compact wavepacket impinging on the QPC from 
the left. Subsequent time evolution then leads to the wavepacket heading through the QPC, with par- 
tial reflection immediately from the region of the QPC if ti{E) is less than unity. The probability of 
immediate reflection, due to the QPC itself, is 

^rnmediate ^ ('^]^\ f | (^) dE (8) 



y m 

where (Ti{E) = J2i ISieiE)]"^ is the reflection probability of the clean QPC from the left, i.e. ai{E) + 
Ti{E) = 1. The time evolved thermal wavepacket is 



M^,y,t) = I dE ariE) g-W/^+Mi^) ^l;,(x,y,E) 



(9) 



Far to the right, ipi{x, y, E) — > tii(E')-(/^°"*(x, y, E), where ip°'^*{x, y, E) is a unit flux outgoing wave 
fai- to the right, beyond all obstructions. Then after a time when the wave has vacated the scattering zone, 
the norm of the wavefunction on the right is easily shown to be 



) I dE \aT{E)\' \tu{Er = ^ldE (^^) n{E) = Gt, (10) 

i.e. the thermally averaged conductance is obtained by propagating a single coherent wavepacket from 
the left and checking how much of it remains on the right. 

An equivalent approach is to integrate the flux returning through the QPC over time, or the flux 
heading to the right. For incoming channel I, from the left, the wavepacket is, on the left, 

i^T{x,y,t) = I dE ariE) e"^^*/^' i^e^^'^^^^Mv) + X] 5,,,,(ii;)e-^'=^'(^)^V^,(y)| (11) 

On the right, it is (in the asymptotic region beyond obstructions etc.) 

Mx,y,t) = J dEariE) e-'^'/^Y.^,y{E)e"'r'iE)x^^,^y^ (^2) 

r' 

Using the flux operator j = h/in lm.{ip^Vip), in either case we easily recover Eq. ^|for the con- 
ductance. The time average of the flux operator using thermal wavepacket thus gives the conductance. 
The flux approach can be extended to calculate the thermally averaged cross conductance matrix ele- 
ments Gpq^T in a multiterminal an^angement by averaging the flux out of terminal q given it was sent in 
terminal p. 
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The approach works because 1) time averaging removes arbitrary phases between states of different 
energy, 2) there is no degeneracy since by assumption only one mode can pass through the QPC and 
emerge to the right, and 3) we can construct a simple wavepacket whose amplitude in the eigenstates is 
exactly 

w{E) = ^-dfT{E)/dE = {AkBT)-^''' Sech[(^ - EF)/2kBT] (13) 

and 4) the amplitude aT{E) is nearly a Gaussian peaked at Ep, (see below) leading to a compact thermal 
wavepacket. It should be clear from the derivation that if there are two open transverse modes in the 
incoming QPC, then we need to run two independent wavepackets, one in the lowest and another in the 
next transverse mode, etc. We may also compute for example the average wavefunction on the right by 
time averaging the wavepacket, etc. 

thermal wavepacket idea is especially well motivated for Fermi gasses with an applied bias voltage 
between leads, where the relevant weighting function w{E) is approximately a Gaussian, peaked at 
the Fermi energy. This means thermal wavepacket is itself nearly Gaussian, which is more localized the 
higher the temperature, and which has a well defined mean velocity and dispersion in velocity. (Although 
the Fourier transform giving the thermal wavepacket is not a closed form, there is no requirement to 
approximate the wavepacket by a Gaussian; the Gaussian approximation is valuable intuitive purposes 
but the exact thermal wavepacket may be found numerically if desired). The width of the wavepacket is 
just the thermal length: It =0.4 microns at 1.7K, and 0.16 at 4.2K in a sample with Ep = 0.016ey and 
Up = 2.86 X 10^ m/s. The thermal length divided by the wavelength is IpjXp = Ep/ir'^kBT ~ 5 in this 
example. The momentum uncertainty of this thermal wavepacket, 6u/up ~ ks T/2Ep is only about 
1% of the average momentum, so that in a 10/im round trip the wavepacket gains only about O.l^um in 
width at 4.2K. Higher temperatures will shorten the thermal wavepacket, and increase its momentum 
uncertainty. The phase coherence length i^f9\ provides an upper limit to the temperature for practical 
use of the thermal wavepacket, since the £^ needs to be longer than ip', this is assured if fc^T << Ep, 
which is normally not a limitation. Journeys longer than iff, will be of course not be coherent. The real 
part of a thermal wavepacket at 4.2K and Ep = O.OlGef in GaAs/AlGaAs is shown in Fig.Q 

3 Physical picture 

The thermal wavepacket approach permits a powerful intuitive picture and a simple estimation tool for 
understanding interference effects and electron choreography in small devices. The basic scenario is as 
follows: Launched from just to the left of the QPC, a spatially nan^ow wavepacket emerges to the right, 
fanning out as part of an annulus, then suffering small angle scattering, splitting up into pieces due to 
collisions with large and small objects, walls, etc. (Fig. |2li. The pieces will have the same width in 
their direction of travel as the original wavepacket (unless the scattering is resonant and therefore time 
delaying). As the pieces arrive at the original QPC (or another terminal) the rule is simple: separate 
pieces of the wavepacket must arrive at the same time if they are to interfere. (If the pieces exit to a lead 
with several modes open, they interfere mode by mode). If an object which the wavepacket encounters 
is moveable, we may follow the oscillation (in conductance) as the interference of the objects' phase 
changes from constructive to destructive and back to constructive as it is moved in certain directions. 

We now apply the thermal wavepacket to explain several situations, including recently published and 
previously unpublished experimental results. 
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Figure 1: The real part and absolute value squared of the thermal wavepacket at 4.2K in a GaAS-AlGaAs 
heterostructure interface. 




Figure 2: thermal wavepacket emerging from the QPC encounters two impurities which are approxi- 
mately the same radius from the QPC; the amplitude from these can interfere at the QPC (right); the 
amplitude scattered from impurities encountered later cannot interfere with the amplitude from the first 
two, since the arrival of that amplitude will be too late. 
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3.1 Adiabatic and nonadiabatic leads 

Consider a wave traveling in a smooth coherent wire, in the lowest transverse mode. Let the profile of 
the wire open up slowly, as in a trumpet bell, forming one half of a QPC. Intuition suggests that nearly 
100% of the wave will emerge without reflection. One argument for this is that the wavelength is short 
compared to the size of the "bell" , and if "slowly" means small changes over a wavelength there is 
hardly any reflection. Calculations bear this out. (It is interesting that sound waves however have no zero 
point restrictions; the wavelengths can be on the order of the length of the trumpet, much larger than the 
bell. This causes significant nonadiabatic reflection at the bell, especially in the lower registers). The 
emerging, spreading, um-eflected quantum wave defines a reflectionless incoming channel tpQ , by time 
reversal. An adiabatic QPC with one open mode has transmission unity (assuming adiabatic bells on 
both ends and a smooth waveguide between them) and G = 2e^ /h. 

Suppose now that one moveable obstruction is encountered some distance past the QPC. We con- 
sider an adiabatic QPC and a perfect ballistic conductor with no backscattering other than the moveable 
obstruction. The wave returning from this obstruction will have some component in the ip^ channel, and 
this will transmit back through the lead without reflection. Other channels will be 100% reflected with 
no transmission. Neglecting for a moment any double or multiple bouncing off the obstruction, no other 
returning flux exists, and thus no interference pattern will result as the obstruction is moved away from 
the QPC, even at very low or temperature. Since the neglect of multiple bouncing is a good approxi- 
mation at long distances, the prediction is that the fringes will disappear at larger radii, even at very low 
temperature. 

3.2 Fringing 

In a previous publication's^, conductance fluctuations were reported in a QPC device as a charged scan- 
ning probe microscope (SPM) tip was scanned across the sample at long range from the QPC. The 
system was free of obstructions on both sides of the QPC (except for unavoidable impurities and donor 
atom potential flctuations) There were thi^ee striking aspects of this data. First was its high spatial res- 
olution, which was explained as due to the "glint" effect from SPM tip divot, much smaller than the 
divot itself II 31 . Second, the flow was found to be strongly branched, a phenomenon which arises from 
a concatenation of caustic structures from the small angle deflections caused by donor atom potential 
fluctuations j3j El- Third, the electron flow branches were decorated with interference fringes, even very 
far from the QPC. These fringes are robust to thermal averaging, which at first is surprising. For exam- 
ple, if the fringing was due to interference from imperfections near the QPC and the backscattering from 
the tip, they should have died out in well under a micron from the QPC, due to averaging over different 
wavelengths in the thermal average. 

We can easily explain the persistent fringing using the thermal wavepacket picture. As the electron 
wavepacket emerges from the QPC, it encounters the small impurities, which each scatter an s-wave 
in all directions as the outgoing wave passes by. (See Fig,|2l). Some of this amplitude is directed back 
towards the QPC. We assume (as is the case in the experiments) that the impurities are sparse enough that 
multiple scattering from the impurities is a negligible contribution at the QPC. The SPM tip also scatters 
the wave coming from the QPC in many directions, although not uniformly, as it is a larger target. Some 
of this amplitude also heads back toward the QPC. (If the SPM "divot" is several wavelengths across, 
then the scattering in a given direction can be attributed to just a part of the surface of the divot- the 
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"glint" effect.) As the SPM amplitude returns to the QPC, any amplitude which is also returning after 
scattering off impurities can interfere. It is thus immediately cleai^ that only those impurities within a half 
a thermal length of the SPM tip are able to effect fringing. The net impurity amplitude arriving at the 
same time as the SPM tip amplitude is a fixed background against which the tip amplitude beats as the 
tip is moved radially away from the QPC. Since the thermal wavepacket is typically many wavelengths 
wide in its direction of travel, the main effect of moving the tip is to cause a phase shift of its amplitude 
and thus regular fringes in the conductance. Larger radial shifts of the SPM (on the order of the thermal 
length) mean that the SPM amplitude is interfering with amplitude from a new set of impurities. This 
picture explains why the fringes seen in the experiment "aim" back to the QPC (i.e. the normal to the 
fringe antinodes points back to the QPC). Constant radius movement of the tip does not change the phase 
of its returning amplitude at the QPC, and thus no fringes are crossed in this direction. According to 
this picture, there should be no die off of the fringing at large radius, for distances much less than i^, 
except for the geometric fall off of the SPM tip signal, which decays as where r is the radius from 
the QPC. However the average signal also dies off at this rate, so the fringes remain unabated in relative 
strength at large r, which is in fact observed. The dynamics of small angle donor potential scattering 
leading to branching leads to interesting corrections to this picture, which can also be explained by the 
thermal wavepacket model. 

Impurities whose main contribution is to backscatter the direct amplitude coming from the QPC must 
therefore reside in an annulus of mean radius R, of width of the order of the ballistic thermal length in 
order to contribute to the interference with the tip at radius R from the QPC. This explains the fringing 
in the first generation of experiments ||3j|, where weak impurity backscattering interferes with the SPM 
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Figure 4: (left) Snapshot of a thermal wavepacket returning from one bounce off the distant mirror, and 
two bounces off the close mirror. A node has developed from the interference between the two, seen as 
a vertical black region above the center of the QPC, This wave has difficulty propagating through the 
QPC, since the excited transverse mode is reflected (close channel), (right) This images is taken at the 
identical moment as the one on the left, with no changes except the outer mirror has been moved about 
1/4 of a wavelength. Now an antinode develops above the QPC, and much more flow back through the 
QPC occurs, corresponding to lowered conductance. 

amplitude to cause fringing. 

Stronger scatterers such as a concave mirror, however, can reflect considerable amplitude back to 
the walls containing the QPC, from which it can bounce back to be reflected again. Such a path, which 
has a length twice that of the single bounce, will only be able to interfere with the direct backscattering 
from the SPM tip if it is positioned approximately twice as far away as the mirror, so its direct reflecting 
path has the same length as the double bounce. At this distance, the single and triple mirror bounce 
interference arrives too early or too late to interfere with the first bounce from the more distant mirror 
(Fig. Ill , and the double bounce fringing emerges essentially uniquely at this distance. This was observed 
in experiments in which the second mirror was the SPM divotl^l2J. 

Figure |4] shows the squai^e of a thermal wavepacket which has bounced off two miiTors, in two 
different cases. In both cases the wavepacket had earlier emerged from the QPC and was then, roughly 
speaking, cut into two pieces by the first mirror. One piece, which hit that mirror, returned toward the 
wall, and the other piece continued its journey toward the second mirror. The first piece, after bouncing 
off the wall and hitting the first mirror again, returns toward the QPC, aiTiving at the same time as the 
first bounce from the second mirror. This is the time at which this figure shows the wavefunctions. They 
interfere with essentially the classic two slit interference pattern, corresponding to their different paths, 
which results in a sequence of nodal lines approximately in the direction of travel. The flux back through 
the QPC is strongly affected by the alignment of the nodal lines with respect to the QPC. The nodal lines 
remain essentially fixed during the arrival of the two packets. If an independent run is made with either 
of the two mirrors slightly shifted, the nodal lines shift. The two cases shown differ by a shift of the outer 
mirror by about 1/4 of the mean deBroglie wavelength. Because the wavepacket interacts twice with the 
close mirror, the rate of shift of the fringes is twice that of the outer mirror. 
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Another scenaiio obtains when the divot is close to the QPC: The amplitude from the first bounce off 
the divot interferes with paths which have bounced twice off the divot; this effect was held responsible 
for the fringing seen in imaging of emerging transverse modes near a QPC14|. The first return bounce 
creates an amplitude against which the double bounce can interfere; since the double bounce aquires a 
phase relative to the single bounce at the same rate the single bounce does to any fixed amplitude, the 
double bounce fringes look "normal", with a fringe spacing of half a Fermi wavelength. However, the 
double bounce scenario is wiped out by thermal averaging if the total pathlength to the obstruction and 
back is greater than a thermal length It, about 0.4 microns in 0]. The thermal wavepacket approach 
shows why: the single bounce amplitude will amve back at the QPC and clear away before the double 
bounce begins to amve if the obstruction is more than a thermal length away. The single and double 
bounces will then be unable to interfere 

Now suppose the QPC has internal "defects", or the energy is such that the electron has to tunnel 
through it (no open channels), making it nonadiabatic. The reflection from the divot into scattering chan- 
nel ijjQ can now travel back thi^ough the lead (with partial reflection), emerging on the left to interfere 
with the amplitude that was reflected from the QPC as it approached from the left side. If the obstruction 
is far enough away that the initially reflected amplitude on the left escapes the QPC region to the left 
before the backscattered amplitude from the obstruction arrives, then once again there can be no inter- 
ference between them and no fringing will result. Suppose this is the case at some temperature T. Then, 
at some lower temperature, fringing will appear- due to this scenario, as the thermal wavepacket becomes 
more extended. 

4 Quantum control with electron phase 

The effective time resolution tt in the thermal wavepacket (the time it takes to pass a given position) 
changes with temperature, is It/^f = h/TrksT. In the familiar case of 4.2 K at Ep = 16 meV, this 
is already half a picosecond; at 20 K, tt = 10~^^s. Events separated by times greater than tt behave 
classically, while quantum interference remains for events within tt of eachother. This permits filtering 
of interference "noise" as distinct from the "signal", by arranging the signal to amve from separate paths 
with Tt- 

A simple realizable experiment exploits this time resolution, illustrates the power of the thermal 
wavepacket viewpoint, and points the way to a new coherent electron switch. Suppose we construct a 
Michelson interferometer involving a 50% splitter mirror and two 100% reflection mirrors in the canon- 
ical arrangement (Fig. [Sjl. This scenario illustrates a coherent control switch for the electron pathway. 
When the outer 100% reflecting mirrors are at nearly equal distance, the contrast is nearly 100% for left 
versus down control; the contrast lessens until when one mirror has been moved a thermal length (divided 
by 2), after which no left-down contrast remains. The key notion is that both wavepackets returning to- 
ward the splitter must amve at the same time if they are to interfere most effectively. Delays in arrival 
of one compared to the other mean that the leading edge of one wavepacket, and the trailing edge of the 
other have no "partner" with which to interfere. Enough movement of one mirror to cause the arrivals at 
the splitter to be nonoverlapping result in no selection for left versus down, and four wavepackets being 
produced. 

Assuming that the mirrors are equidistant, the Michelson scenario points to the possibility of a 
"warm" but still phase controlled, coherent, electron switch. A gate lying above one of the two paths 
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Figure 5: A schematic of a feasible Michelson interferometer experiment with electrons. The thermal 
wavepacket is shown impinging on the 50% miiTor in a. In b, the split packets are heading for their 
respective mirrors, and in c, they are about to interfere coherently after experiencing equal length round 
trips to give d, the constructive interference for the "left" as opposed to "down" paths. In e-h, the right 
mirror has been moved out, so that the wavepackets return at different times and are unable to interfere, 
spliting into four pieces with no "left-down" selection. The thermal length (divided by 2) determines the 
minimum distance the right mirror needs to move to remove the interference. 
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taken by the electron wavepacket after scattering off the splitter could phase shift one of the split packets 
relative to the other, using a small charge and with almost no change in amval time. This would change 
the electron path from down to left, for example, after the packets recombine at the splitter. 

5 Simulating higher temperatures 

To demonstrate the validity of the thermal wavepacket picture that has been presented, we have imaged 
the electron flow from a QPC using a low-temperature SPM in a configuration similar to the two mirror 
geometry discussed previously. We lithographically defined a metal gate that acts as one of the mirrors, 
and our movable tip acts as the other. Figure 6a shows a schematic of the device geometry and the 
area of the scan, and subsequent panels show the tip scan data at increasing "temperature" (see below). 
Interference fringes appear in the image from distinct round trip trajectories of the same length, to within 
a thermal length. Increasing the electron temperature decreases the thermal length, and localizes the 
extent of the fringes, Fig.Q. Note the corresponding size of the thermal wavepacket, which shrinks with 
increasing temperature. 

To increase the electron "temperature", we used an arbitrary function generator to generate a time 
varying voltage signal across the QPC. The voltage from the function generator injects electrons of that 
energy through the QPC, and mimics the energy distribution of electrons at the intended temperature. 
The amount of time at any given voltage is given by the probability of the electron having that energy at 
the desired temperature. Therefore, the voltage spends the most amount of time at V, which corresponds 
to the Fermi energy and the smallest amount of time at the highest and lowest voltages, which are the 
tails of the distribution. As the temperature is increased the voltage spends more time away from V 
because there is a lai^ger spread in energy of the electrons. 

6 Conclusion 

Thermal electron interferometry has strong analogies to while light interferometry. The simplest ex- 
amples are two slit interference and weak localization coherent backscattering. The central peak in the 
two slit case, and the weak localization phenomenon are examples of interference phenomena which 
survive averaging over different wavelengths. We have introduced a new way of calculating and under- 
standing thermally averaged electron interference in semiconductor heterostructures. A single "thermal 
wavepacket" whose width in the direction of travel is the thermal length can yield thermal averages of 
measurable quantities such as conductance. Various experiments were described in terms of thermal 
wavepacket dynamics, including some earlier data, and some reported here for the first time. A proposal 
for a coherent Michelson switch was also given. 

Computing thermal wavepackets can be an advantageous way to obtain thermal averages numerically. 
The alternative is to obtain fixed energy solutions and average over them explicitly. Although we do not 
anticipate here what methods might be employed to make the fixed energy calculations efficient, the 
efficiency of fast Fourier transform propagation for wavepacket dynamics is legion ll4ll . Obtaining a 
theoretical SPM image, which includes all the aspects implied by the presence of the divot, is a time 
consuming task, requiring a complete scattering calculation and thermal average for every new position 
of the SPM tip. This gives one pixel of information for each tip position. Even a modest image of 1000 
pixels requires 1000 independent, thermally averaged runs, the thermal wavepacket approach enjoys 



12 




Figure 6: (a) Schematic diagram showing the location of the images of electron flow. Images of electron 
flow at three temperatures, (b) 1.7 K, (c) 4.2 K and (d) 8.4 K showing the robustness of the fringes near 
the radius of the arc. See the text for an explanation of the procedure used to simulate these temperatues. 
The white line in (b) shows the location of the one-dimensional line sample used in Figure 
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Figure 7: Graph of the measured signal as a function of the distance from the arc for a series of simulated 
temperatures, along the strip shown in Fig.|6l The strength of the signal decays more quickly away from 
the arc for increased temperatures but remains robust at the arc position. In each case, a black outline 
shows the dimension of the corresponding thermal wavepacket 

many advantages in this situation, one of which is re-use of the calculations up to the time when thermal 
wavepacket hits the divot. 

Acknowlegement: We thank Jacob J. Kritch for a critical reading of this manuscript. This work 
was supported at Harvard University by DARPA grant DAAD19-01-1- 0659 and by the Nanoscale Sci- 
ence and Engineering Center (NSEC) under NSF grant PHY-01 17795; and NSF CHE-0073544; work 
at UC Santa Barbai^a was supported by the Institute for Quantum Engineering, Science and Technology 
(iQUEST). 

References 

[1] B. J. van Wees et al., Phys. Rev. Lett. 60, 848 (1988). 

[2] D. A. Wharam et al., J. Phys. C Sohd State Phys. 21, 209 (1988). 

[3] M.A. Topinka et al.. Nature 410, 183 (2001). 

[4] M.A. Topinka et al.. Science 289, 2323 (2000). 

[5] J. A. Katine, et al., Phys. Rev. Lett. 79, 4806 (1997). 

[6] M. Blaauboer and E. J. Heller "Statistical distribution of Coulomb blockade peak heights in adia- 
batically pumped quantum dots" Physical Review B 64241301/1 (2001). 

[7] R. Landauer, IBM J. Res. Dev. 1, 223 (1957); M. Bttiker, Phys. Rev. Lett. 57, 1761 (1986). 



14 



[8] S. Datta, "Electronic Transport in Mesoscopic Systems", Cambridge, 1995. 

[9] G.F. Guiliani and J.J. Quinn, Phys. Rev. B 26, 4421 (1982). 

[10] M. A. Topinka, R. M. Westervelt, and E. J. Heller, "Imaging Coherent Electron Flow", Physics 
Today, p. 47, December 2004. 

[11] S.E.J. Shaw, R. Fleischmann, E.J. Heller, lcond-mat/0 105 3541 (2001): S.E.J. Shaw, Ph,D, Thesis, 
Harvard University, 2002. 

[12] B.J. LeRoy, A.C. Bleszynski, K.E. Aidala, R.M. Westervelt, A. Kalben, E. J. Heller, S. E. J. Shaw, 
K.D. Maranowski and A.C. Gossard, submitted. 

[13] M.A. Topinka, R.M. Westervelt, and E.J.Heller, "Imaging Electron Flow", Physics Today, Dec. 
2003. 

[14] R. Kosloff,in Numerical Grid Methods and Their Application to Schrdinger's Equation, edited by 
C. Cerjan (Kluwer Academic, Dordrecht, 1993), p. 175. 



15 



